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The energetics of virtually all binding function* in proteins 
is the culmination of a set of molecular interactions. For 
example, removal of a single molecular contact by a point 
mutation causes relatively small reductions (typically 0.5-5 
kxal/mo!) in the free energy of transition-slate stabilization 
(for reviews see Fersht (1987) and Wells and Estell (1988)]. 
protein-protein interactions (Latkowski et ah, 1983, 1989; 
Ackers & Smith, 1985), or protein stability [for review see 
Matthews (1987)] compared to the overall free energy asso- 
ciated with these functional properties (usually 5-20 kcal/mol). 
Thus, it is possible to modulate protein function by mutation 
at many contact sites. In fact, to design large changes in 
function will often require mutation of more than one func- 
tional residue. 

There is now a Urge data base for free energy changes that 
result when single mutants are combined. A review of these 
data snows that, in the majority of cases, the sum of the free 
energy changes derived from the single mutations is nearly 
equal to the free energy change measured in the multiple 



jimple additivity breaks down. The first is where the mutated 
residues interact with each other, by direct contact or Indirectly 
through electrostatic interactions or structural perturbations, 
so that they no longer behave independently. The second is 
where the mutation causes a change in mechanism or rate- 
limiting step of the reaction. It is important to note that the 
additive effects discussed here do not change the molecularity 
of their respective reactions. When the molecularity of the 
reaction changes [as in comparing the free energy of binding 
of one linked substrate (A-B) versus the sum of two fragments 
(A plus B)]. large deviations from simple additivity can result 
from entropic effects (Jencks, 1981). Although the focus here 
is on enzyme activity, similar conclusions may be drawn from 
mutations affecting protein-protein interactions, protein-DNA 
recognition, or protein liability. Some practical examples and 
applications are discussed. 
ADornvrrv Relationships 

The change in free energy of a functional property caused 
by a mutation at site X is typically expressed relative to that 



of the wid-type protein as AAGpQ. Such free energy changes 
for two tingle mutants (X and Y) can be related to those of 
a double mutant (designated X,Y)byeq 1 (Carteret aL, 1984; 
Ackers & Smith, 1985). The AG, term (also called the 
AAGfjcy, = AAGpq + AAGry, + AG, (1) 
coupling energy. Carter et si, 1984) should reflect the extent 
to which the change in energy of interaction between sites X 
and Y affects the functional property measured. It is possible 
for AG, to be either positive or negative depending upon 
whether the interactions between the mutant side chains reduce 
or enhance the functional property measured. Furthermore, 
the AG, term should not exceed the free energy of interaction 
between side chains at sites X and Y except in caseswhere 

first applied to evaluating the functional independence of 
residues mutated in tyrosyl-tRNA synthetase (Carter et al., 
1984). In one case the sum of the AAG value* for tingle 
is equal to thai of the double mutant, indicating 



a large discrepancy, suggesting the rites were interacting. 

Simple ADDrnvm in Transition-State Bimxno 
Interactions 

The strength* of noncovalent interactions are strongly de- 
pendent upon the nature of the two groups and- the distance 
\r) between them. For example, the free energy of charge- 
charge, random charge-dipole, random dipole-dlpoic, vai 
Waats attraction, and repulsion decay as 1/r, l/r\ !//*, 



and I/r", respectively [for review see Fersht (1985)]. 1 
when the side chains at sites X and Y are remote to one 
another and assuming no large structural perturbations, the 
AG, term should be negligible and eq 1 thus simplifies to 
AAGkv, as AAG m -I- AAG m (2) 
This situation, here referred to as simple additivity, Is generally 
observed except where side chains arc dose to each other or 
when one or both of the mutants change the rate-limiting step 
or reaction mechanism. These principles are well illustrated 
from data of additive mutational effects on tranaition-sute 
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another mutant, R198M (Scruiton ct al, 1990). These basic 
residues are not in direct contact, but both side chains form 
a alt bridge with the ^-phosphate group of NADPH. Indeed. 



D»u .retaken from Table I and repre«e n t rrnttanUi hunt subtilisin 

fivould not contact one another. The dashed fc« . *>« rf l. 
and the solid line is a best fit to all the data. 



in trai 



nergy (AaG T ) 



caused by a mutation can be calculated from eq 3 (Wilkinson 
et al., 1983), in which * is the gas constant, T is the absolute 
(W*ttUu- 

temperature, *„, is the turnover number, and K v is the Mi- 
cfcaeHs constant for the mutant and wild-type enzyme against 
a fixed substrate. AACf represents the change in free energy 
to reach the transition-state complex (ES») from the free 
enzyme and substrate (E + S). 

To analyze the proposition that the interaction energy terra, 
AC™, is relatively small when the sites of mutation (X and 
Y) are remote to one another, AAGj values were collected 
from the literature where side-chain substitutions in the 
multiple mutant are beyond van der Waab contact (>4 A 
distant) from each other (Table I). There are at least 25 
examples distributed across five different enzymes where 
MGi values can be calculated for the individual and multiple 
mutants assayed in at least two different ways. Among these 
are examples where electrostatic interactions, hydrogen 
bonding, and steric and hydrophobic effects have been altered 
separately or in combination with others. The X-ray structures 
of the wild-type proteins show that the wild-type side chains 
are not in contact. Modeling suggests the mutant side chains 
are beyond possible van der Waals contact unless the mutant 
side chains were to cause significant changes in the overall 
protein structure. Such large changes are rarely observed in 
structures of site-specific mutant proteins (Katz ft Kossiakoff. 
1986; Alber et aL, 1987{ Howell et aL, 1986: Wilde et aL, 
1988) or even highly variant natural proteins (Chothia ft Leak, 
1986). 

A collective plot of the sun^ the M^vahics forte 

mNerT gi^Trein«rtably strong correlation (** - 0.92) 
wth a Senear unity (Figure 1). Theshnpfattaerpretatto 
is that the interaction term, AGV, is small compared to the 
overall effects on AAGfovr It is formally possible that there 
are Urge and compensating effects between side chains X and 
Y that systematically lead to small net values for AfiV 
There are some notable exceptions that weaken the corre- 
lation within the dau^Tabk^. ^ g^^^redurtaw 
gives a less than additive effect, especially when combined with 



with NADPH as ct 

oftheAA<7 T vi ... 
mutants in subtilisin (D99K and E156K) overestimates the 
effect of the multiple mutant when assayed with sn'Arg but 
not with a Phe substrate (Russell ft Fersht, 1987). Such, 
discrepancies are not too surprising because charge-charge 
interactions fall off as I/i- and can exhibit long-range effects 
in proteins [for example, see Russell and Fersht (1988)]. The 
physical basis for other large discrepancies not involving 
electrostatic substitutions is less dear but may involve unex- 
pectedly large structural changes or changes in enzyme 
mechanism (see below). 

These additivity tests are not particularly dorninated by one 
of the single mutants in the sum. The average contribution 
(±SE) for the most dominant mutant in ech sum calculated 
from the 69 additivity tests given in Table I is only 68% 
(±15%) of the total sum (theoretical is -50%). Furthermore, 
the plot in Figure 1 is not analogous to graphs of correlated 
variables, where A is plotted versus the sum of A + B, because 
in Figure 1 the values cut the yaxis are determined Inde- 
' i those on the x-axis. 



Complex AoDmvnv in Transttion-Statb 
Stabilization— When AGfa * o 

(A) Change In Interaction Energy between Situ XandY. 
Where residues Xand Y are dose enough to contact, ft is nwre 
likely that the AGfa term will be significant, There are 1 1 
examples collectively from tyrosyl-tRNA synthetase and 
subtilisin that fit this category (Table II). 

A series of mutants in tyrosyl-tRNA synthetase at positions 
48 and 51 (Carter et al., 1984; Lowe et aL, 198S) show com- 
plex additivity (Table II). His48 and ThrSl in the wild-type 
structure are next to each other on adjacent turns of an er-henx. 
His48 hydrogen bonds to the ribose ring oxygen of AIT while 
ThrSl can make van der Waals contact with ATP. TheTSlP 
mutation increases the catalytic efficiency of the enzyme in 
some assays by more than -2 kcal/rao! (Wilkinson et al, 
1984). However, when this mutation is combined with mu- 
tations at position 48, the effects are not simply additive. An 
X-ray structure of the T51P mutant indicates there are no 
structural changes in the a-heiix (Brown et al„ 1987). Instead, 
it is suggested that the TS1P mutant b improved over wild 
type because the wild-type enzyme contains a bound water in 
the vicinity of Thr51 that disfavors substrate binding. Blow 
and co-workers (Brown et aL, 1987) argue that the change 
in solvent structure propagated to position 48 may account for 
the complex additivity. In the previous section, the double 
mutant (H480.T51A) exhibited nearly simple additivity 
(Table I). Presumably, the smaller and less hydrophobic 
alanine substitution at position 5 1 should not introduce as huge 
a change in solvent structure as the rKnro&ione ring cf proline. 
In the case of subtilisin (Table II). Glnl56 is near toe top 
of the PI binding crevice while Glyl66 is at the bottom. In 
the wild-type enzyme these sites do not make direct van der 
Waals contact, but large side chains substituted at position 
IWcanlwmoc^tocwitact the residue at position 156. In 
fact. X-ray structural analysis shows that an Asn side chain 
at position 166 makes a good hydrogen bond with Gtal56 
(Bott et aL, 1987). Moreover, all of the substitutions are polar 
or charged, the energetics of which are rapected to be the most 
long range. Thus, the mutant xide chains alter substantially 
■ ■ L - een positions 156 and 166. 
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Tible I: Comptmon of Sums of AAGr* fwm Component Mutxrts v. tht Multiple Mut«t 
AitfContiaOaeAnotkcr 


Where the Ms 


UM or Wild-Type Side 


Chains A> 






AAGr' 




cedent mutants- .... 




WUltt Ml 





TyiosytoRNA Synthetase 



ATP/PP, 
ATP/tRNA 
Tyr/PP, 
Tyr/lRNA 



MTX 



C35Q + H4SG* 






+1.20 +1.04 




+130 


-i-tjQS +1.13 


+118 


+1.68 


+1.14 +1.12 


+126 


+132 


+0.31 +1.12 


+U5 


+1.20 


C33G + TJIP 






+1.20 -1.91 


-0.71 


-1.14 


+1.05 -2J3 


-1J0 


-1.8* 


+1.14 -0.64 


+O50 


-0.7* 




+082 


+0.21 


C3SG + TS1C» 






+IjOS -0.93 


+0.12 


-0.22 


+1.14 -0.91 


+023 


-0.13 


H48N + T51A* 






+0.26 -0J8 


-012 


+0.04 


-013 -032 


-04S 


-037 


T40A + H45G' 






+5.02 +3.1S 


+8.17 


+6.95 


+S.13 +2.44 


+7.57 


+6.67 


Ret Trypsin 






G2I6A + G226A' 






+2.75 +3.13 


+5.88 


+5.07 


+2.19 +4.91 


+7.10 


+5.90 


DihydrcMate Reductase U 
F31V + LS4C/ 






+1.6 +19 • 


+4J 


+4.3 


♦12 +19 


+5.1 


+4.5 



E156S + 
-1.43 
-060 
-0.15 
+1.70 
-086 
-0.61 
-024 
EIS6S + 
-1.43 
-O60 
-0.13 
+ 1.70 
-086 
-061 
-024 
EIS6S, 
Y217L* 
-1.67 
-0.96 
-0.53 
+1.33 



Y2I7L + 0169A' 

-087 -0.62 

-0J6 -0.32 

-041 -027 

-0.08 -0.30 

-OJ2 -039 

-0.29 -0.«6 

-012 -0.41 
Y217L 
-0.87 
-0.36 
-041 
-0.08 
-0.32 
-0.29 
-0.12 



-2.30 
-0.96 
-056 
+1.62 



-062 
-0.32 
-0.27 



6+B156S' 
7 +0.T7 
-0.62 



+087 
-1.41 
-1.17 
-0J9 



+1.33 
-l.U 
-0.S4 
-0.32 



-1.14 
-0»2 
+047 
-1.41 



D99K + E156K 
+1.29 +2.12 
+0.13 -0.49 



0169A.+ SSJ? 
Y2I7L 

-1.46 +4.96 
-1J33 +4.40 
S24C. 

»^ + 0166A 
Y2I7L 

+4 J I -O40 
+196 +094 
S24C. E1S6S, 
H64A,+ G!69A. 
G166A Y2I7L 
+4.11 -1.46 

S24C, + GI66A. 

H64A G169A. 

Y217L 
+4.96 -1.76 
+4.40 +0*2 

E. eo*f Gluuthkme 
A179G + RI98M* 

-1.10 -062 
+O08 +168 

AI790 + R2041 

-1.10 +0.41 
+O08 +142 

R198M + R204L 

-062 +041 
+168 +142 

a nor. + 
A,79Q + R204L 

-1.10 -05I 
+O08 +3,70 
RI98M + 



-062 -1.54 

+168 +047 
■U04L + A179G, 
IUWL+ R198M 

+0.41 -1J2 

+142 +111 

R179G + RI98M + R204L 

-1.10 -0.62 +041 

+O08 +168 +142 



... iu of ATP-depemfan pyrophosphate ewhwge (ATP/PPj) or iRNA tiuugjng (ATP/tRNA) 

under ulariting coikBUom for tyrosine end rke vera for Tyr/PPj «ckra|e end Tyr/OWA chargtoe, ♦ et el (1 985). ^^™fpP£& 



•Carter etal. (1984). The assays refer » mi 



aetel (1986). 



vetta assay refers wfcewticei of tyrosyl . 
The ATP/Tyr ud Tyr/Tyr ectivetion eneyr refer la formation of tyrwyl 
* J /J^r« lyrc«iiie« B J AlT.rei ) »«ively. 'Crete; et si. (1985). The sebttrits we. i>-V,|.Uu-(X)h 
• "erLys(K)orAxg(R>. 'Meyer eteL (1956). The ligand «i either cShydrc/oktr (H,F) or mrtho 

'he Pi (X) residue (Schccbter A . 

ht (1987). The mhurstc was benzoj+t-Vel-Gly-L-Arg-s-nitroejiillde 
). 'Carter etal (1989). The wtatnle wis i«ccinyl-L-Pbe-L-AU-l.-HU-t-(X)-e- " 



at tuccinyl-L'AJl*L- 

Ah (A), Lye (K). Met (M). Pfc* (F), or Tyr 0 

«lctwyl-L-AuM.-AtaH-Prc-l-r^ 



rin where the Pi reside* (X) 

te(MTX). 'Webs etti (1987a). The 



either Pfc (F) or Tyr (Y). 'Scretton et U. (1990). The eray foltowd the reduction of oridncd liuuthjorie by WAPH or NADPH. 
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T.bk II: Conparison of Sums of iA&r* from Cwpsaeat 14 
vi die Multiple Mutsat Where the M«unt Side Chains Can < 



ATP/PP, 
ATP/tRNA 

Tyr/lRNA 
ATP/Tyi* 
Tyr/ATP 

ATP/Tyr 



ATP/lRNA 

ATP/Tyr 
Tyr/Tyr 
ATP/lRNA 

ATP/Tyr 
Tyr/Tyr 
ATP/lRNA 



TyrosyHRNA Synthetase 
H48G + T5IF* 
+1.04 -1.91 
+1.13 -2.35 
+1.W -064 
+1.12 +O50 
+0.9S -1.99 
+1X7 -038 
H4IN + T5IP 
+0.18 -19* 
+OW -0.38 
-O02 -2.23 
N4S0 + TJ1P 
+037 -0.94 
+0.41 -1.00 
+1.26 -1.05 
Q4SG + T5IP 
-1.09 



-047 
-1 .11 
+048 
+1.63 



-0.57 
-0.59 
+0.21 



-2X5 



-1X5 



-1J7 
H48Q + T51P 
+226 -1.99 
+3.13 -038 
+311 -2.23 

Subulum BPN' 
EI56Q + GI66D' 
-1X4 +1.27 
-045 +1.83 
+2.15 +043 
E156S + GI660 
-0.39 +1.27 
-085 +1.83 
+1.6* +053 
E156Q+G166N 
-1.71 -0.11 
-1X4 +0.14 
-045 +018 
+115 +C48 
EI56S+GI66N 
-1.44 -0.11 
-0.59 +0.14 
-OSS +0.18 
+1.68 . +0.48 
EI56S + G166K 



+027 
+2.75 
+0.8S 



+0-23 
+1.38 
+2.68 

+048 
+OM 

+122 



-0.90 
-0.27 
+2.73 



-0.4J 
-0.67 
+2.16 



+1X7 
+0.77 
+1X2 
+0.17 
+1.04 
+0.82 



+086 

+0.45 
+0.90 



+1.17 
+1.48 
+1.26 



+075 
+0.16 
+0.26 

+0.74 
+0.66 
+049 



-0.51 
-0.85 
-0.78 
+1.26 



r 

r 





SSa.;* 




> 










+ -4 -2 0 2 4 



now 1 DataintatafimTii&IIiariiwiuitiafnbfili^M} 
of unity dope tad the wild lint rtpraeaa the bat fit. 



the component (ingle mutants (** - 0.72). Moreover, the 
slope of the line (0.61) is much below unity. This indicates 
that the function of one residue is <»mpromised by mutation 
of another. Of tie 40 additivity example*, the average con- 
tribution of the most dominant single mutant totheitimof 
the values is 71% (±13%) of the total. Thus (as in 
figure 1), both single mutaau can contribute tubstantiaUy 
to free energy changes measured in the multiple mutant. 
However, this data set is derived from mutations at only two 



In summary, complex additivity can be observed when 
mutations at sites X and Y change the intramolecular inter- 
action energy between sites. Tbiscan be mediated by direct 
steric, electrostatic, hydrogen-bonding, or hydrophobic in- 
teractions or indirectly through large structural changes in the 
protein, solvent shell, or electrostatic interactions. Complex 
additivity is most likely to occur where the sites of mutation 
are very close together and larger or chemically divergent aide 



-1.44 



-1.03 



-OS» 
-083 
+1.68 +041 
E156Q + G166K 



-1.71 



-1.03 



-2X7 
-1J2 
+2X6 



-095 
-1.12 
+1X8 



M -0.45 -IJ7 
K +2.15 +0.51 t«j» 

•See Tabk I for desen^Uoa assays. * Lowe el iL (1985). 'Carteret 
»L(I984). 'Weill el al (1987b). 

In these six examples there are large and systematic discrep- 
ancies between the sum of the AACf values for the tingle 
mutants and those of the corresponding double mutant (Wells 
etaL, 1987b). In abncstall cases, the sum of the AACjvaiues 
for the single mutants it much greater than the value for the 
multiple mutant. Nonetheless, the AAG} value predicted from 
the sum of the single mutants does have the same sign as that 
for the double mutant, so that the single mutants predict 
qualitatively the effect on the multiple mutant 
A plot (Figure 2) of the collective data set from Table II 
st to that seen in Figure 1. The AAG T values for 
re poorly with the sum of 



(*) Mutations as Sites X or Y Change ike Enzyme 
Mechanism or Rale-UmltlngSlep. If the catalytic functions 
of two or more residues are interdependent, then a mutation 
of one residue can affect the functioning of the other(s). This 
form of complex additivtty is well illustrated form 
the catalytic triad and oxyanion bindin " 
(Carter ft Wells, 1988, 1990). In the o . 
of sobtilisin (Figure 3), the rate-limiting step in amide bond 
hydrolysis ii trawfer of the proton from Ser221 to His64 with 
nucleopailic attack upon the Kissile carboayi carbon. This 
is accompanied by electrostatic stabilization of the protonaied 
imidazole by Asp32 and hydrogen tending to the oxyamon 
by the side chain of AsnlSS and the main-chain amide of 
Ser221. Mutational analysis shows that once the catalytic 
Ser221 is mutated to Ala (S221A), additional mutations in 
the triad or oxyamon binding site cause no further loss in 
catalytic efficiency (Table III). 

The 8221 A enzyme retains a catalytic activity that ts still 
10* above the solution hydrolysis rate (Carta ft Wells, 1988). 
It is proposed that this residual activity is derived from re- 
maining transition-state binding contacts outside of the cat- 
alytic triad coupled with solvent attack upon the earbonyl 
carbon from the face opposite position 221 (Carter ft Wells, 
1990). This proposal is based on a. model showing that there 
is no room for a water molecule near Ala221 once the nmstrate 
is bound. Furthermore, conversion of Asnl55 to Gly enhances 
the activity of the S221A mutant by -1.2 kcal/mot (Table HT). 
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AJM55/ 




Tibhlll: , ... - 

wthc AACt* far Multiple Mataaia la the Catalytic ' 
Oxytaion Banm, Site af SubUlMa — 



S22IA + KMA* 
+8.93 +144 

saiA + DJZA 
+s\»3 +fc52 

H64A + D32A 
+8JM +4J2 

S12tA+H64A + D3IA 
+193 +M4 +M2 +H.29 



+11.78 +«.*> 
+I5.4J +IW 
+tSJ» +7*8 



+8JI. +M3 + 
SniA + NlSSO' 

+1.93 +3.01 + 

'All auymcs were mjuyed with the jubrtmtc wceiayl-L-AU.L-AU- 
L-Pro-t-Phe-pHrint-niMi. »Ctrtcr >ad WtSl (!««)• 'Carter tod 



+WS 
+7.70 



This Is consistent with the opposite-face solvent attack 
mecnaiusmorS#IA,l>ec»iise theoxyaiuon (Figure 3> would 
develop away from Asnl55 and the NI55G mutation improve* 
solvent accessibUUy to the sdstttc carbonyl carbon. 

Complex additiviiy is also seen for subtilisin mutated at 
positions 64 and 32, The double (H«4A4»2A) and corre- 
sponding single mounts show a linear dependence upon hy- 
droxide ion concentration (between pH 8 and 10) that may 
reflect hydroxide assistance in the deprotonatton of the Or 
of Ser22l (Carter ft Wells, 1988). Thus, once His64 is 
converted to Ala; AspM b a liability, presumably by elec- 
trostatic repulsion of hydroxide ion. rNotetnc-1.3kcal/iwi 
improvement in AACJ for the double mutant (H64A.D32A) 
compared to H«4A alone, Table III.] 

In summary, if an enzyme mechanism relics upon cooper- 



mutations within this subset can result in large values for 
AG^i). In (act, if the mechanism is changed substantially, 
residues that were a catalytic asset can become a liability. 
Simple additivity can also break down when-one or more of 
die nutations cause a change ia t^rate-limiting step. In an 

that enhance the activity, but the cm 



activity could not go beyond 
(Albery ft Knowfcs, 1976). 
Additive Effects on Substrate Binding 
The analysis above considered changes in Wnding free cn- 



b(E + S) to yield 
Thest 



the bound transitica w m t e c omp lex (&S*). 1 
kinetic analysis for subtiUsin and tyrosyMRN Aiym^tetase is 
such that the K u values approximate the enzyme-aubstrate 
dissociation constant.*;. AiLm^aiialysi! fcanwdemeum* 
latioos ofAAC^B,, (from K u values) or MO m (ten** 
values) yields qualitatively the same result* (not shown) as 
sbowin Tahks landU a«l Fajnres 1 and 2. Tbtis, deviations 



the energetic to form the &S complex or those to reach BS». 
Aoomvs Effects on fwrnat-naim htxnMxma 

Tbe fust dear examples of additive binding effects caused 
by amino acid replacements in proteins were reported by 
Laskowski et aL (1983) and reviewed by others (Ackers.* 
Smith, 1985; Horovkz ft Righi, 1985). C ' " 



have been isolated and sequenced from the eggs of different 
bird specie* (Empie ft LaskowsU, 1982; UskowsU et at, 
1987). This is a nested set of proteins because for any one 
0/ them avian uuWduw 

one or a few amino acid substitutions. Moreover, thcassd- 
datksscoojnu«fJCJcfto 

prtteittasmvaiycwanencnBomrang^ Lasfowskj 
etaL (1983, 1989) have shown that the effect of a |>en residue 
replacement on X. is about the tame irrespective of the in- 
hibitor scaffold the replacement is made in. 
In addition 10 ovomucoid, four additivky examples h»ve been 



tetrameric hemoglobin (Ackers ft Smith. 1985). Three ad- 
ditivity examples have been analyzed for interactions of hOH 
with its receptor (B. C Cunningham and J. A. Wells, un- 
published results) and one example for association of synthetic 
variants of the iWase S peptide with RNase S protein 
(Mitchittsoo ft Baldwin. 1986). The entirety of due data set 
is not tabulated because much on the ovomucoid inhibitors 
and bOH is unpublished. Nonetheless these reaearchers were 
kind enough to provide their data formatted so it could be 
plotted collectively in Figure 4. These dau consist of 91 
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FKXAE4: Plot thowing the sum ofchanees in free energies of toadies 
at prote'm-protcin intcrfacct for component mutants versus the 
corresponding multiple mutant. Data represent bteractiana between 
ovomucoid third omnia end various serine proteases (a) (R.Wynn 
and M. Uskewski, personal communication), regulatory htwrsce 
ofoAbemoglobis (•) (Ackers & Smith. 1983). hGH and W receptor 
(Kippted A) (B. Cunaisghim and J. Wells, posonal comtmiaicat&n), 
andRNase S peptide and S prateia (■) (Mitcbiwon A Baldwin. 
19«6). The daihed line represent* a line of unity slope, and the solid 
line is the best fit. 

kcal/md). The plot shows a very strong linear correlation (& 
m 0.96) with a slope near unity. Although the data for the 
at sorted to evaluate changes at ' 



andallofthe other examples represent nonoonuct sites. Thus, 
the large date base derived from natural variants of ovomucoid 
third domain, as well as a smaller number of examples from 
several other proteins, indicates that multiple mutations at 
protein-protein interfaces commonly produce simple additive 
effects. 

Additive Effects in DNA-Protein Interactions 

One of the clear advantages in analyzing DNA-protetn 
interactions is the abiSty to apply powerful selections that make 
analysis by random mutational studies feasible. Additivity 
in DNA-protetn interactions was first demonstrated by re- 
version analysis of X repressor (Nelson & Saner, 1985). A 
mutation that decreased the binding affinity for the X operator 
site (K4Q) was reverted by mutations at several second sites 
(E34K,G48S,andE83K). When these second-site rrverunti 
were introduced Into wild-type X repressor, they caused in- 
creases in affinity similar to those observed in the first-site 
suppressor mutant (K4Q). 

Functional independence for mutations at DNA-protetn 
contacts has been demonstrated by additive effects for mutants 
of CAP (catabolite gene activator protein) and its operator 
sequence (Ebright et ml., 1987) as well as toe repressor and 
its coaesponding operator sequence (Ebright, 1986). Simple 
additivity of mutational effects in the operator sequences for 
Cro repressor (Takeda et al., 1 989) and X repressor (Sarai ft 
Tekeda, 1989) has been most systematically demonstrated. 
Simple additivity has also been reported for multiple mutations 
in the toe repressor (Lehmlng et al., 1990). In tact, simple 
additivity is so predictable in DNA-protein interactions that 
the observation of complex additivity has been used to predict 
specific DNA-protetn contacts in the foe repressor -operator 
complex (Ebright, 1986). 
Additive Effects on Protein Stability 

The first systematic analysis of additive effects of sitc- 
spodifo mutations on protein stability was reported by Sbortle 
and Meeker (1986). Five multiple mutants in staphylococcal 



VWL + G79S + GMV 



N-Terauaal Domain of X Repressor 
G46A + G4IA* 
It +0.7 +0.9 +145 



13C + C54V 



+1.2 



-0.7 



IJC + CJ4T 
+U +03 
I3C + CS4T + R96H 
t +1.2 +0.3 -2J 
I3C.C54T + R96H 
al melt +U -24 

I3C + CS4T + A146T 
+IJ +0.3 -1.5 
13C.C54T + AI46T 



KringleJoftPA 
H64Y + WW 
+19 +0-7 



•Sbortle and Meeker (1984). »Heet« et at. (1916). 'Weteet el al 
(I9IS). 'Sandbars and TerwiBlger (1919). *R. KeUey. penwaal 
cemimialcatioi. 'OtSfwikl aad Uskowskl (1990). M5-CHO refen 
loiglycoiyUiionofA^i. »Herle et at (I9»0. 

I from a group of random single 
i initially for their ability to affect 
the stability of the enzyme in vivo. The component mutants 
t with each other in the muttipk 



also been determined for the N-ter m i n al d om a i n of X repressor 
(one example: Hecbt et al„ 1986), the o-si»burnt of £.co/f Trp 
synthetase (one example; Huifc et aL, 1986), T4 lysozyme (mix 
examples; Wetzel et al., 1988), the gene V product of bac 
' e fl (one example: Sandberg ft Terwilriger, 1989), 
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combined, is one of the most powerful tooii in drjugning 
functional propcrtk* in proteins, m approach baa been 



Plot showing mm of change* is fret aver© of anfo 
r em mutants and resulting multiplt mutant. Data are 
uw, .able IV ind represent ruphykxxccal nuclease (■), N-termin&l 

synfteUsep). mda^Bnerep^.toc^lineoraraty 
slope, and the solid fine represent! the belt fit. 

natural variants of ovomucoid third domain (two examples; 
Otlcwstd ft Laskowsld. 1990), and the KrIngle-2 domain of 
human tissue plasminogen activator (t-PA) (one example: R. 
Kelley, personal communication). 

Collectively, this data set gives a high linear correlation {ft 1 
• 0.94) and dope near unity (Figure 5). The generally iimpfc 
additive behavior is somewhat surprising given the highly 
cooperative nature of protein folding. Then are discrepancies 
in some of the activity examples besides the staphylococcal 
nuclease mutant (V66L,G88V). For example, the U 
kcal/tnol discrepancy for the Y175C.0271E double mutant 
in Trp synthetase (Table IV) is proposed to result from the 
fact that these residues are in direct contact (Hurle et al., 
1986). Furthermore, proximity effects may account for the 
large differences between the sum of the component mutants 
and the multiple mutants for the o-hefical double glycine 
mutant G46A.Q48A in X repressor (Hechl et al.. 1986). and 
when combining R96H with the CJ-C97 disulfide mutant in 
T4 lysoxyme (Wetzel et al., 1988). In contrast, an exchange 
or two side chains that contact one another (V35I and I47V) 
in the hydrophobic core of the gene V product of fl phage 
produced simple additive effects (Sandberg ft TerwilKger, 
1989: Table IV). It should be noted that tins data base ex. 
htbiting simple additivity may be biased for single mutants 
that stably fold, because severely unstable proteins are more 
difficult to express. 

By analogy to transition-stale binding effects, one can 
certainly imagine Instances where the stabilizing effects of 



at high temperatures can become controlled by a chemical step 
such as deamidation (Ahem et al., 1987), so that additional 
mutants that stabilize the folded form of the protein may be 
irrelevant. Another obvious example where panmiex additivily 
can be observed in protein stability fa the stabilizing effect of 
disulfide bonds and itoitcovalent intramolecular contacts that 
require interactions between two or more residues. In these 
cases, the stabilizing interaction between two side chains can 
be broken with only one mutation. 

Applications of Additivity in Rational Protein 



mactiration, such as X repressor (Hecht et aUm wfcflfam 
(Bryan et al, 1987; Cunningham ft Wells, 1987; Pamohano 
et aL, 1989), lunamycin nucleotidyltransferase (Liao et at, 
1986; Matsumura, 1986). neutral protease (Imanaka et aL, 

1986) , and T4 tysozyme (Wetzel et al.. 1988; Matsumara et 
aL, 1989). This strategy has been applied to enhancing the 
catalytic efficiency of a weakly active variant of subtffisra 
(Carter et aU 1989), engineering the substrate specificity or 
tubtihsin (Wells et aL, 1987a.b; Russell ft Fersht, 1987) and 
the coenzyme specificity of glutathione reductase (Scrutton 

et aL, 1990), designing protease inhibitors with exquisite • 
protease specificity (Laskowsld et aL. 1989). and recruiting 
human prolactin to bind to thehGH receptor (Cunningham 
et at, 1990). In addition. additMty principles have been used 
to engineer the pH profile of uibdlfain (Russell ft Fersht, 

1987) and to design the affinity and specificity of X repressor 
(Nelson ft Saner, 1985). 

For this approach to work does not require that ail the 
component mutants act in a simply additive manner but just 
that their effects accomuiate. For example, despite the oom- 
ptaaddWvhy of effects in the catalytic triad of aibtiUsin, there 
are mutagenic pathways that are energetically cumulative for 
installing the triad (Carter ft Wells, 1988; Wells et aL. 1987c). 
Starting with the triple mutant S221A.H64A.D32A, there is 
a progressive enhancement for installing Ser221 (-1.1 kcal/ 
mot), then His64 (-1.0 kcal/mol), and finally Asp32 (-6.5 
koal/mol). Another cumulative pathway of Ser221, then 
Asp32. and finally Htt64 is possible if the Ser22l,Asp32 in- 
termediate were to use HisP2 substrates (Carter ft Wells, 
1987). Elaborating such cumulative pathways is important 
for understanding how a catalytic apparatus may have evolved 
and is practically useful for considering how to install «ch 
catalytic machinery into weakly active catalytic antibodies. 

Conclusions 

In the majority of cases, combination of mutations that 
affect substrate or transition-state binding, protein-protein 
interactions, DNA-protein recognition, or protein stability 
exhibits simple additivity. Simple additivity fa commonly 
observed for distant mutations at rigid molecular interfaces 
such as in protein-protein and DNA-protein interactions, 
where the mutations are unlikely to alter grossly the structure 
or mode of binding. 

Large deviations from simple additivity can occur when the 
sites of mutations strongly interact with one another (by 
making direct contact or indirectly through electrostatic in- 
teractions or huge structural perturbations) and/or when both 
sites function cooperatively (as for the catalytic triad and 
oxyanion binding site or luhtilisin). Changes at sites that can 
contact each other do not always lead to complex additivity; 
this may reflect relatively weak interactions between the two 
sites or indicate that the interactions are compensatory and 
appear to be weak. 

It is important to point out the magnitude of errors in 
predicting the free energy effect m the multiple mutant from 
the component single mutants. Generally, for those eases 
exhibiting simple additivity (Figures 1, 4, and 5), the dis- 
crepancy ia free energy between the sums of the components 
and multiple mutant* is about ±25%. Part of this is the result 



the real 

Nonetheless, this means that if the total free energy change 
is about 3 kcal/mol, the change in the equilibrium constant 
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(related by - 10-V"" - 155) wiH often be off by a 

factof of 4. Thus, while the free energy effects accumulate, 
significant deviations will occur in predicting the final equi- 
librium constants when component mutants contribute a large 
free energy term. 

Simple additivity reflects the modularity of component 
amino adds in protein function. This results from the fact 
that the perturbations in energetics and Structure resulting 
from most mutations are highly localized, to the past six years, 
live mutagenesis strategy has been extremely effective 
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abstract: Femtosecond spectroscopy was used in combination with site-directed mutagenesis to study the 
influence of tyrosine M210 (YM210) on the primary electron transfer in the reaction center of Rhodobacter 
sphaeroides. The exchange of YM210 to phenylalanine caused the time constant of primary electron transfer 
to increase from 3.5 ± 0.4 ps to 16 ± 6 ps while the exchange to leucine increased the time constant even 
more to 22 ± 8 ps. The results suggest that tyrosine M210 is important for the fast rate of the primary 
n transfer. 



The primary photochemical event during photosynthesis of 
bicteritthtorophyll- (Bchl-) containing organisms is a light- 
induced charge separation within a transmembrane protein 
complex called the reaction center (RC). The crystal struc- 
tures of RC's from Rkodopstudomonas (ftps.) oMdls and 
Rhodobacter (Rb.) sphaeroides have been solved to high 
resolution (reviewed in Ddsenhofer and Michel (1989). Chang 
et aU (1986). Tiede et al (1988). and Rees et «L (1989)]. The 
RC from Rb. sphaeroides contains three protein sabunits 
referred to as L, M, and H. according to their respective 
mobilities in SDS-polyacrylamide gels. Associated with the 
Land M subunits are the cefaclors, consisting of four Bchl 
a, two bacteriopbeophytin (Boh) a, one atom of non-hemc 
ferrous iron, two quinoaes (Q A and Qg), and in some species 
one carotenoid [reviewed in Parson (1987) and Feher et aL 



(1989)]. The cefaclors are arranged in two branches (Figure 
I) with an approximate C; axis of symmetry. The kinetic data 
support a model in which the primary electron transfer pro- 
ceeds after light absorption by the primary donor (a special 
pair of Bchl referred to as P, reviewed in KimuuerandHolten 
(1987)]. The absorption of light generates the excited elec- 
tronic state P*. which has a lifetime of approximately 3 ps. 
An electron is transferred from P along only one branch (the 
so-called A-branch). It is generally accepted that after ap- 
proximately 3 ps the electron arrives at the Bph on the A-side 
(HJ and after 220 ps it reaches Qa. The role of the accessory 
Bchl located between P and H A (referred to as has not 
been definitely assigned. Recently, we have shown that at 
room tonperature an additional kinetic (r - 0.9 pa) component 
is detectable (Holzapfel et al„ 1989). The spectral properties 
and the kinetic constants lead to the conclusion that the 
corresponding intermediate is the radical pair P*B A " (Hol- 
zapfel et al.. 1990). 
Additional Intriguing points concerning the process of 
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